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A  frequently  observed  ship  wake  feature  in  Synthetic  Aperture 
Radar  ( SAR )  images  is  a  dark,  narrow  line  along  the  ship  track. 
These  features  have  been  tamed  turbulent  wakes,  although  it  is  net 
clear  that  turbulence  alone  is  responsible  for  their  appearance,  it 
has  been  suggested  that  vortices  proceed  by  the  ship's  hull  may  be 
responsible  for  the  suppression  of  surface  waves  near  the  sh’p  t>*ac< 
ano  an  enhancement  of  the  waves  near  trie  edge  of  the  s.-ectned  area. 


The  primary  purpose  of  '  the "  study  reported  here  was  to  use  a 
hydrodynamic  model  for  ship-generated  vortices  as  input  to  existing 
radar  backscatter  mode’s,  and  to  use  this  hybrid  model  to  produce 
simulated  SAR  ship  wake  signatures  for  comparison  with  actual  mea¬ 
surements.  For  the  limited  number  of  cases  considered,  the  simulated 
wake  signatures  agreed  closely  to  those  observed  in  actual  SAR  data. 
Qur  results  also  indicated  that  the  SAR  signatures  are  related  to  tne 
vertex  Currents  under  lew  wind  conditions  and  to  the  Strain  nates 
associated  with  the  snip-generated  vertices  under  high  winds. 


Additional  studies  inc'uded  the  racar  resolution 
SAR  snip  wa<e  features  as  well  as  examination  of  the 
form  of  a  snip  wake  image  for  information  which  could 
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INTRODUCTION 

Ship  wake  patterns  exhibiting  a  wide  diversity  of  appearance  have 
been  noted  on  Synthetic  Aperture  Radar  (SAR)  images  collected  with 
various  systems  under  a  broad  range  or  env ironmental  conditions.  As 
a  result  of  a  careful  examination  of  the  available  data  and  a  cata¬ 
loging  of  these  observations,  in  addition  to  several  recent  theo¬ 
retical  analyses,  it  appears  that  the  wake  structures  within  the 
existing  oody  of  observational  data  can  be  classified  into  tnree 
general  categories  of  phenomena  as  shewn  in  figure  i.  ’’hese  cate¬ 
gories  include:  ;1)  surface  waves  generated  directly  by  the  snip, 
(2)  turbulent  wakes  or  vortices  containing  relatively  persistent  but 
non-propagating  currents,  which  are  visible  by  SAR  through  the  inter¬ 
action  of  ambient  waves  with  these  Surface  currents,  and  (3)  inter¬ 
nal  waves  generated  by  the  ship,  which  again  are  visible  because  of 
their  interactions  with  short  surface  waves.  ERIM  is  actively 
involved  in  all  three  phases  of  $AR  ship  wake  studies.  An  earlier 
study  (Lyden,  et  a!.,  1985a),  carefully  examined  the  SAR  imaging 
characteristics  of  sh ip-generated  surface  waves .  That  study  ; r. pl¬ 
eated  tnat  the  brignt,  narrow  7-wakes  deserved  behind  ships  ur.cer 

low  wind  conditions  could  be  adequately  explained  by  a  f irst-erde”* 
Sragg  scattering  mechanism.  ERIM  has  also  performed  research  into 
Ship-generated  internal  waves  (Lyden,  1985b). 

"he  purpose  of  the  study  reported  here  has  been  to  investigate 
the  narrow  dark  bands  frequently  observed  in  SAR  imagery  extents* ng 
berrnd  moving  ships.  These  features  occur  when  the  ocean  sjrface  is 
sufficiently  rough  to  yield  a  measurable  radar  return,  wnicm  '-s 
suppressed  in  tne  region  near  the  snip  track.  The  suppress  *  cr.  *s 

frequently  more  pronounced  at  .-band  (23.5  cm  wavelength)  tnan  at 

.(—ban  a  (3.2  cm  wavelength) ,  arc  is  sometimes  accompanied  by  a  tr; gn: 
line  on  one  or  pot.n  edges  of  the  dark  area.  An  example  of  t.nese 
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Figure  i.  Overview  of  Shin  wake  Generation  Mechanisms 
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features  is  shown  in  Figure  2  for  simultaneously  obtained  L-  and 
X-band  imagery. 


The  features  described  above  have  been  termed  turbulent  wakes 
(Shuchman,  et  al.,  1983),  although  it  is  not  clear  that  turbulence 
alone  is  responsible  for  the  appearance  of  these  features.  The 
effects  of  turbulent  currents  on  surface  waves  are  not  well  under¬ 
stood  and  the  persistence  of  the  phenomenon  is  somewhat  problematic. 
It  has  been  suggested  that  horizontal  vortices  produced  at  the  edges 
of  the  ship's  hull,  as  well  as  a  net  rearward  velocity  within  the 
turbulent  wake,  may  be  responsible  for  the  suppression  of  surface 
waves  near  the  ship  track  and  an  ennancement  of  the  waves  near  the 
edges  of  the  smoothed  area  (Swanson,  1984),  These  characteristics 
are  shown  graphically  in  Figure  3.  The  vortex  mechanism  would  also 
help  explain  the  persistence  Of  turbulent  wake  features  in  SAR 
imagery.  Turbulent  wake  features  have  been  observed  extending  up  to 
20  km  aft  of  the  generating  ship  in  Seasat  imagery. 

In  this  study,  a  first  attempt  is  made  at  using  a  hydrodynamic 
model  *or  these  vortices  as  input  to  existing  radar  backscatter 
models,  and  to  use  this  hybrid  model  to  produce  simulated  SAR  wake 
signatures  for  comparison  with  actual  measurements.  Presented  below 
is  a  description  of  the  hydrodynamic  and  radar  backscatter  models  we 
employed.  This  is  followed  by  a  detailed  quantitative  comparison  for 
a  single  data  set  which  appears  to  validate  the  application  of  this 
model . 

The  model  was  then  exercised  to  predict  the  expected  SAR  signa¬ 
ture  for  several  hypothetical  ship,  SAR,  and  environmental  condi¬ 
tions.  These  results  can  be  used  to  aid  in  the  design  of  an  optimum 
SAR  system  for  the  detection  of  these  wakes.  This  is  followed  by  a 
presentation  of  several  SAR  images  which  appear  to  qualitatively 
agree  with  model  predictions.  Also  presented  in  this  report  is  an 
investigation  into  the  effects  of  SAR  resolution  on  ship  wake 
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X-Band 


Figure  2.  Simultaneously  Obtained  L-  and  X-Band  Optical 
SAR  Imagery  of  a  Turbulent  Wake  Feature 
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SURFACE  CURRENTS  DUE  *0  VORTICES 
(Typically  i  *  i GO  cn/si 


NET  rearward  VELOCITY 
(TYPICALLY  1  -  10  cn/s) 


■<L'LL  VORTICES 
NOT  ”C  SCALE ) 


Figure  3.  Diagram  of  Turbulent  Wake  Structure  Showing  Possible  Role 
of  Vortices  and  Net  Rearward  Velocity  Within  Turbulent 
Region  (after  Swanson,  1984) 
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detection  as  well  as  a  study  on  the  use  of  the  Fourier  transform  of  a 
$AR  ship  wake  image  for  detection  applications.  Finally,  a  conclu¬ 
sions  and  recommendations  section  is  presented.  Included  as  an 
appendix  is  a  discussion  on  the  detectability  of  surface  patterns  in 
$AR  imagery, 
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SAR  MODELING  OF  VORTEX  WAKES 


In  this  section,  a  description  of  the  hydrodynamic  model  for  the 
vortex  wakes  is  given  as  well  as  a  description  of  the  radar  back- 
scatter  model  we  utilized  in  our  simulations.  This  is  followed  by  a 
quantitative  comparison  of  a  simulated  wake  signature  with  actual 
measurements  for  a  case  where  the  various  ship  and  environmental 
parameters  were  well  known.  A  series  of  simulations  were  also  per¬ 
formed  to  help  determine  an  optimum  SAR  configuration  for  the  detec¬ 
tion  of  vortex  wake  features.  Several  SAR  images  of  turbulent  waxe 
features  are  also  presented  which  qualitatively  support  tne  vertex 
mechanism.  Finally,  a  brief  . ivestigation  of  the  effects  of  SAR 
resolution  on  ship  wake  detection  is  presented,  as  well  as  a  discus¬ 
sion  on  using  the  Fourier  transform  of  a  ship  wake  image  for  detec¬ 
tion  purposes. 

2.1  HYDRODYNAMIC  MODELS 

The  hydrodynamic  aspects  of  our  study  can  be  grouped  into  two 
categories.  The  *irst  is  the  manifestation  of  surface  currents 
through  vertex  generation.  The  second  is  the  interaction  of  tnese 
currents  with  the  snort-scale  3ragg  waves  which  produce  the  r-jca^ 
backscatter.  Each  of  these  topics  is  discussed  below. 


2.1.1  VORTEX  WAKE  MODEL 


'he  components  of  the  proposed  vertex  wax e  mcce1  are  i  ’  lustr  ats 


■n  iaure  A. 


,n  this 


igure,  a  ship  is  shewn  travel 


ng 


irreuan 


water  at  steacy  s-ceec  u. . 


awa  v 


)rwaro  momentum. 


It  '•adiates  a  Kelvin  waxe.  which  carries 
'he  turbulent  waxe  produced  by  tne  $no 


therefore  must  have  a  net  rearward  momentum,  as  snewn  in  tr.e  f 'ours . 
In  aeditiun,  a  vortex  pair  is  produced  by  the  snip  null,  whicn  causes 
a  chance  in  toe  1 ateral  flew  field  behind  the  ship.  And  final iv, 
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when  a  strong  thermoc 1  i ne  is  present,  the  flow  field  of  the  vortex 
pair  can  interact  with  it,  causing  pressure  and  velocity  changes  at 
depths  which  generate  internal  waves.  All  of  these  hydrodynamic 
phenomena  give  rise  to  horizontal  currents  at  the  ocean  surface  which 
interact  with  surface  waves  to  produce  changes  in  surface  roughness 
which  are  potentially  detectable  by  $A5.  The  present  report  con¬ 
siders  only  the  effects  of  the  vertex  pair. 

A  ship  generates  a  flow  field  due  to  a  vortex  pair  which  is  pro¬ 
duced  and  shed  by  the  hull.  In  all  displacement  ships,  the  water 
•nust  accelerate  as  it  passes  arourd  the  bottom  of  the  snip  hull,  as 

illustrated  in  Figure  5.  This  gives  rise  to  a  pressure  drop  along 
the  bottom  of  most  of  the  hull,  due  to  the  Bernoulli  effect,  this 
effect  has  been  confirmed  on  ship  hulls  by  numerous  measurements,  and 
gives  rise  to  a  negative  lift,  or  downward  force  which  is  well-known 
to  naval  architects  as  "sinkage-in-tr im"  or  "squat."  This  force 
causes  tne  ship  to  rice  lower  in  the  water  at  higher  speeds. 

This  force  L  can  be  described  approximately  by 

I  (  M  A  0  r  Z  /  l  \ 

w  =  \^o,  j  ~pvs  , 


wh  ers 


is  the  ship  beam,  i  its  length, 


the  water  density,  -j  the 


snip  speed,  and  C  a  coeffici ent  between  0 
P 

hull  snape.  The  coefficient  C  has  been  estimated  as 

p 


iecencrrg  on 


a  1.5; 


\  3/2 


U/ 


:  ? ' 


Crcm  the  theory  of  lift,  it  can  be  shewn  that  tnis  force  must  give 
r;se  to  the  snedding  of  at  least  one  vortex  pair,  and  the  total  cir¬ 
culation  r  in  tne  vertex  must  be 


odor 


VV 

7 7~ 


w  r.ere  f  is  a  constant  approximate : y  equal  to  2. 
distribution  pf  li*t  on  the  null  pf  the  Ship. 


arising  *rcr  fe 
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Pressure  Drop  Due  to 
Bernoulli  Effect 


Figure  5.  Schematic  Showing  How  Velocity  Increases  and  Pressure 
Decreases  8eneath  Hull  of  Displacement  Ship  (from 
Swanson,  1984) 
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The  definition  of  circulation  r  Is  that  It  is  the  total  inte¬ 
grated  vorticity  per  unit  length  along  track  in  one  of  the  vortices. 
When  the  vorticity  in  such  a  vortex  is  highly  localized  in  the  vortex 
core,  the  velocity  flow  field  due  to  the  vortex  car.  be  described  very 
simply  as  a  tangential  velocity  field  V&  by 


♦*» 


where  r  is  the  circulation  in  the  vortex,  and  r  is  the  radius  to  the 
vortex  center. 


the  vortices  shed  by  the  hull  will  have  the  geometry  and  orienta¬ 
tion  shown  in  figure  6.  This  will  produce  upwasn  in  the  center  of 
the  wake  behind  the  ship,  and  generate  a  spreading  or  diverging  flew 
in  the  wake  center  at  the  surface.  This  lateral  surface  velocity  '$ 
designated  V  . 

Because  the  boundary  condition  for  the  water  surface  is  that  the 
*low  must  be  primarily  horizontal,  a  pair  of  image  vortices  are 
required  above  the  surface  to  satisfy  this  condition.  The  inter¬ 
action  of  the  null  vortices  with  their  images  causes  a  lateral  motion 
of  the  null  vortices  away  trem  one  another,  "he  latera'  spreading 
speed  7„  of  the  hull  vortex  core  will  be 


here,  z  is  the  equilibrium  depth  of  toe  vortex.  It  can  be  shown  that 
eacn  vortex  comes  to  an  equilibrium  depth  of  about  d/2,  which  is  very 
close  to  the  depth  at  which  > t  is  produced,  and  consequently  its 
motion  after  shedding  is  primarily  horizontal. 


"he  res^ 


1 ater al  velocity 


■  e  i  i 


at 


he  s j r r a c ; 


given  by 
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where  c  is  the  lateral  separation  distance  between  the  vortices,  Z  is 
the  depth  of  the  vortices,  and  y  is  the  lateral  distance  from  wake 
centerline  or.  the  surface.  As  a  result  of  this  analysis,  we  see  that 
the  vortex  pair  generates  a  lateral  spreading  velocity  on  the  surface 
behind  the  ship,  and  this  current  will  also  have  an  effect  on  the 
amplitude  of  any  surface  waves  which  are  present,  thereby  offering 
possibility  of  detection  by  SAR.  The  interaction  of  these  /ortsx- 
induced  surface  currents  with  small-scale  surface  waves  will  be  dis¬ 
cussed  in  the  next  section. 


2.1.2  WAVE— current 


The  microwave  reflectivity  of  the  ocean  surface  is  primarily  a 
function  of  the  small-scale  (centimeters  and  decimeters)  surface 
roughness.  For  incidence  angles  between  about  20*  and  70*,  Brace 
scattering  is  believed  to  be  the  dominant  scattering  mechanism,  "he 
Bragg  scattering  mode'  yields  a  radar  cross  section  per  unit  area 
(for  the  backscatter  case  with  the  plane  of  incidence  perpend i co 1 ar 
to  the  y-axi s )  of 


4  4 

3,(9',  3  Atk  COS  9: 


s;n  e,  O' 


whe-e  k  is  toe  electromagnetic  wavenumber, 


9  is  the  incidence  angle, 


W ( <  ,  k  )  is  toe  wave  height  spectral 

x  y 

scattering  coefficients  (Valenzuela,  1973) 
tne  backscattered  power  is  proportional  to 
density  at  the  Sragg  wavenumber 


density,  arc  g.^lel  are 
.  Thus,  to  first  Order, 
tne  surface  wave  soectra' 


* 

■<3  =  (2*  s-n  a,  0)  .  -S' 


."or  tre  purposes  of  our  wa<e  stupes, 
of  spectral  density  of  Sragg  waves 
produced  our^-ant  oatta-'n  and  whether 
a  SAR. 


we  are  interested  in  the  change 
at  each  ooint  in  :r,e  vc'tev- 
these  changes  are  ebse-v as ’ a  by 
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To  achieve  this,  a  compute1”  program  has  been  written  and  imple¬ 
mented  on  a  Zenith  Model  100  microcomputer .  This  program  was  ini¬ 
tially  written  for  use  in  internal  wave  studies  and  is  described  in 
detail  by  Lyzenga  (1984).  In  summary,  this  program  solves  for  an 
initial  wavelength  outside  the  current  pattern  corresponding  to  the 
Bragg  waves  at  each  point  in  the  current  pattern.  The  program  cal¬ 
culates  the  change  in  spectral  density  for  these  waves  due  to  current 
interactions.  The  basis  for  these  calculations  is  the  conservation 
of  wave  action  soectral  density  (Phillips,  1977)  which  is  defined  as 


N 


where  F  is  the  wave  heidht  spectrum  and  u  is  the  radian  wave  fre- 

1  yp 

quency  (w  *  ( gk )  j .  The  wave  height  spectral  density  at  equi¬ 
librium  is  assumed  to  be  given  by  the  Phillips  spectrum 

F0(k)«AiT4  (1C) 

so  the  wave  action  spectral  density  at  equilibrium  is 


The  ceviation  *rrn  the  ecuil iorium  state  's  soeci f  *ed  by  tne  scectr a’ 
perturbation  f,  which  is  defined  as 


f 


Nf  k  )  -  N  ('<) 

rpy- 


I 


2  ' 


Note  that  this  perturbation  function  applies  to  the  wave  he’ght  soec- 
trjm  as  well,  thus,  the  perturbation  is  directly  related  to  tne 
expected  change  in  radar  oacxscatter  based  on  a  Bragg  scatte^-rc 
model  -”Eq.  *7;].  To  relate  this  perturbation  value  to  the  mere  com¬ 
monly  used  bac<scattar  change  in  decibels,  one  would  simply  use  the 


re! ationsh  ip 
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Several  points  should  be  clarified  regarding  this  approach,  The 
initial  wave  conditions  are  based  on  a  Phillips  spectrum,  which 
assumes  that  the  waves  are  saturated.  This  may  not  be  true  at  low 
wind  speeds.  Also,  it  is  not  clear  how  much  waves  can  grew  beyond 
saturation  before  the  onset  of  breaking.  This  would  minimize  the 
effects  of  positive  perturbations  calculated  by  cur  model.  It  should 
be  mentioned  that  the  wind  growth/re  1 axati on  rates  used  in  cur  model 
were  based  on  the  work  of  hughes  ('9781.  Finally,  we  should  note 
that  only  transverse  currents  associated  with  tne  ship-generated 
vortices  nave  been  considered.  We  nave  not  included  the  effects  of 
tne  rearward  axial  current  in  this  work. 


2.2  MODEL  RESULTS 

Included  in  this  section  is  a  detailed  comparison  between  our 
model  predictions  and  actual  measurements  made  cn  data  collected 
curing  the  1983  Georgia  Strait  Experiment.  Also  included  are  a 
series  of  model  calculations  of  wake  detectab i 1 i ty  for  a  variety  of 
ship,  radar,  and  env  ironmenta’  parameters.  These  results  can  be  used 
to  determine  an  optimum  sensor  configuration  for  vortex  ...axe  detec¬ 
tion.  Finally,  several  SAP  images  of  ship  wakes  are  presented  arc 
discussed  within  the  context  cf  tne  vortex  wake  model. 

2.2.1  DETAILED  COMPARISON 

Although  there  exists  a  great  deal  of  SAS  imagery  containing 
snips  and  turbulent  wake-like  features ,  there  are  very  few  cases 
where  tne  wake-generating  s.n-.o  parameters  ana  wind  conditions  we^e 
well  documented.  One  exception  is  the  SAR  data  collected  during  tne 
1383  Georgia  Strait  Experiment.  The  required  snip  parameters  for  cur 
model  are  length,  beam,  and  speed.  Surface  w,:nd  speed  and  direction 
are  also  necessary  vnputS .  During  portions  of  the  Georgia  Stra-t 
Experiment,  the  Navy  tug  Cuaoaw  was  making  controlled  ■'•uns  as  cart 
of  a  snip  wake  study.  During  these  •••uns ,  the  ship  speed  was  -e'd 
constant  and  the  surface  wind  conditions  were  monitored. 
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^resented  in  Figure  7  is  an  L-band  image  of  the  Quacaw  and  its 
wake  wnich  was  collected  during  tne  1983  Georgia  Strait  Experiment. 
The  stern  waves  and  turbulent  wake  generated  by  the  Qua  paw  are 
clearly  visible.  The  approximate  wind  conditions  during  this  run  are 
shown  graphically  on  the  image.  O'*  kilometer  aft  of  the  ship's 
location  along  the  line  A— A 1  indicated  on  Figure  7,  radar  backseatter 
values  were  extracted  and  are  shown  plotted  in  Figure  8.  These  mea¬ 
surements  indicate  a  peak  in  radar  intensity  at  the  upwind  edge  of 
the  wake  feature  which  is  about  11  dB  above  the  minimum  intensity 
level  at  the  wake  center.  A  simulation  of  this  wake  was  performed 
correspond ing  to  this  location.  The  parameters  we  used  in  this  simu¬ 
lation  were  consistent  with  tv  Ouapaw ' s  during  the  run  and  were 
based  on  the  fast  vortex  decay  model  discussed  by  Swanson  (1984).  we 
allowed  the  separation  of  the  vortices  [b  in  Eq.  (6)]  to  vary  to 
match  the  width  of  tne  observed  wake  feature.  rhe  surface  current 
profile  for  this  simulation  is  Shown  in  Figure  9(a).  The  correspond¬ 
ing  spectral  perturbations  which  were  simulated  for  this  case  are 
presented  in  Figure  9(b).  These  were  generated  for  a  surface  wave¬ 
length  of  3G  cm  which  is  approximately  the  3ragg  wavelength  for  fe 
data  studied.  Figure  9(b)  can  be  interpreted  as  the  wine  ana  waves 
bcV  traveling  from  right  to  left  entering  tne  current  patter-- 
plotted  in  Figure  9(a).  These  results  indicate  a  strong  positive 
perturbation  at  the  location  of  tne  maximum  adverse  current,  anc  a 
maximum  negative  perturbation  where  the  waves  and  current  are  in  tne 
same  direction  and  the  current  is  i  maximum.  Downwind  cf  this  point, 
the  waves  gradually  return  to  equilibrium. 

A  comparison  between  tne  predicted  and  actual  raoar  intensities 
for  this  wake  feature  is  shown  in  Figure  10.  The  simulated  '•aoar 
intensifies  were  obtained  by  converting  the  spectral  per  turbati  :ns 
shown  plotted  -n  Figure  9(a)  to  equivalent  backseatter  changes 
decibels  using  the  rel  atiensh  ip  given  previously  as  Ed.  (13).  "•'e 

vertical  and  horizontal  position  of  the  :«o  plots  was  adjusted  to 
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Figure  7.  L-Band  Digital  SAR  Irrage  Collected  During  the  1983 

Georgia  Strait  Experiment;  the  Location  Where  Backscatter 
Measurements  were  Performed  is  Indicated  by  the  Line  A-A ' 
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(a)  Surface  Currents  Calculated  for  Line  A-A '  Shown  in 
Figure  7,  and  (b)  Spectral  Perturba ti ons  of  30  cn  Waves 
for  the  Current  Profile  in  (a) 
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Figure  10.  Graphical  Comcari son  of  Predicted  Sadar  Backscatter 
Variation  Based  on  Vortex  Mode1  (dashed)  and  Actual 
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achieve  tne  best  agreement  based  or.  visual  observations.  ?he  close 
agreement  between  the  predicted  and  actual  data  is  obvious. 

Some  notes  of  caution  should  be  made  regarding  these  resuits. 
The  parameters  we  used  1h  cur  mode!  were  adjusted  to  mater,  the 
observed  data,  however,  they  were  very  close  to  what  we  would  have 
predicted  for  the  Quapaw.  The  point  we  are  trying  to  make  in  this 
discussion  is  that  the  vortex  wake  model  discussed  in  the  previous 
section  can  be  used  to  help  explain  ship-generated  waxe  features, 
we  use  the  resu's t  discussed  above  as  justification  for  js'ng  this 
model  to  predict  expected  >*adar  signatures  for  different  ship  arc 
radar  parameters  under  a  variety  of  wind  conditions,  ''his  is  tie 
topic  of  the  next  section. 


2.2.2  MODEL  5-1EOICTIONS 

One  goal  of  developing  the  model  described  above  was  to  be  able 
to  predict  the  detectability  of  vortex  wakes  from  various  ships  for  a 
variety  of  different  wavelength  radars  and  surface  wire  conditions. 
Predictions  of  this  type  would  also  help  define  an  optimum  radar  sys¬ 
tem  for  vortex  wake  detection. 


Two  Ship  cases  were  used  i.n  our  simulations  corresponding  to 
Case  1  and  Case  3  ;n  Swanson  1334} .  *hese  correspond  to  a  smai  1 
f i ve  siloton  displacement  military  snip  and  a  SCO  kilctor  displace¬ 
ment  supertanker .  The  parameters  for  each  of  these  two  snips  are 
given  in  Table  '  . 

"he  decay  characteristics  of  vertex-generated  surface  cu^ssts 
witn  distance  aft  of  the  generating  snip  were  calculated  by  Swanson 
i '  9  S  4  }■  -p-  the  two  snips  we  modelled  and  are  shewn  plotted  -r 
F i cure  1 1 .  It  i s  ’ n teres t * ng  to  note  from  this  plot  that  tne  r53r_ 
ward  axial  current,  a It.nougn  sra''er  initially,  decays  -ere  slow' y 
trar.  the  .'or tex-rgl  ated  lateral  currents  and  eventual  ;y  *'s  ’ arcer  in 
magm  tjde .  Inclusion  of  axia1  currents  in  a  more  comprehensive 
mode ’ 1 ; rq  effort  is  one  recomrenoation  r?Sui t i nq  from  tn  is  work  . 
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SUMMARY  OF  SHI? 


table 

PARAMETERS  ; 


SED  :n  vortex  simulations 


5  <i Teton  Sh * p 
(Case  1 ) 

500  Kiloton 
(Case  3) 

Veloci ty 

TO  m/s 

7.5  m/S 

Length 

TOO  m 

A  00  m 

Beam 

TO  m 

68  m 

Draft 

5  ,m 

26  m 

Displacement 

5  Kten 

500  Ktor 
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Surface  Ships 


"ABAS  DlS/iSiftkj 


500  K  Ton 


5 °0  KTon 

-<f/' 

KTonX 


Voloclty" 


5  K  Ton 


Axlt»  / 

Velocity 


5  *  Ion  Sh‘P  500~K 

CW  D  I 


I  US  •  !0  M/S** 

I  L  *  ’W  M 

j  b  •  1  o  m  j 
*  *  5  M 

Co  ■  07 

•  D"0  *5  K  Ton  : 


K  ton  Ship 
Cm*  3) 

7  5  M/$*c 
400  M 

$8  M  i 

!6  V  I 

07 

500  <  Ton  i 


. . . .  ,oJ 

“JC*y  °f  “ake',"'lucM  S^ace  Currents  as  a  r 

3, stance  9en,nd  Senerat-nn  Sr-  ,  Ct'Cn 

1984)  at‘n9  Ship  (from  s 
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One  set  of  s imu 1 ations  were  generated  for  the  five  kilcton  ship 
corr espording  to  a  location  one  kilometer  behind  the  ship.  The  peak 
latera’  velocity  from  figure  11  is  about  3  cm/s.  We  adjusted  the 
vorticity  in  our  model  to  produce  similar  currents.  The  remaining 
parameters  in  our  model  were  fixed  at  the  values  predicted  by  theory, 
The  surface  current  profile  for  this  case  is  shown  plotted  in 
figure  12.  The  detectab i! ity  of  these  features  as  a  'unction  of 
radar  wavelength  are  shown  plotted  in  figure  13  for  wind  speeds  of 
3,  1.5,  3,  6,  and  12  n/s.  The  surface  wavelengths  we  considered 
incljded  3.2,  5.7,  25,  and  ICC  cm.  At  30*  incidence  angle,  these 

would  correspond  to  the  3ragg  wavelengths  at  X-band,  C-band,  '_-oanc 
and  ICO  cm,  respectively. 

The  parameter  we  have  plotted  in  Figure  13  is  defined  as  the 
absolute  value  of  the  maximum  negative  spectral  perturbation  pre¬ 

dicted  by  our  model.  This  parameter  was  selected  for  several  rea¬ 
sons.  examination  of  historic  SAR  data  sets  clearly  demonstrated 

that  the  turbulent  wake-like  features  we  feel  are  related  to  a  vortex 
mechanism  are  visible  due  to  a  decrease  in  '*acjar  return  (negative 
spectral  perturbation)  aft  of  a  snip  3rd  not  due  to  an  increase 
nadar  backscatter  along  the  edges  of  the  waxe.  Ai>c,  0 m  ~ccei 

assumes  a  saturated  Phillips  spectrum  as  t.ne  initial  unperturbed 
conditions,  it  is  jnclear  new  to  is  spectrum  can  be  enhanced  before 
tne  onset  of  breaking.  The  detectab 1 1 i ty  of  spectral  perturbations 
in  SAR  imagery  is  discussed  in  detail  in  Appendix  A.  For  cur  pur¬ 

poses,  a  spectral  perturbation  of  about  3.15  appears  to  be  tne  thres- 
nolc  of  detectability.  This  value  is  indicated  on  figure  13  by  a 
cashed  line 

"he  results  plotted  in  figure  13  shew  that  two  effects  are  occur- 
-vrg.  when  relaxation  effects  are  not  included  (U  *  C  m/ si,  tne 
sr.ortest  waves  s ; . cw  tne  maximum  pert uoaticn  0 y  the  vertex-gener 2tSd 
surface  currents,  when  relaxation  e'fects  are  included,  the  pertur- 


oat;on  pf  the  snorter  wavelengths  's  sharply  reduced.  ’his  ;s 
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Figure  12.  Modelled  Vortex-Produced  Surface  Currents  One 
Kilometer  Behind  a  5000  ton  Ship  (Case  1) 


radar  division 


Because  waves  of  these  lengths  quickly  relax  back  to  their  equi- 
librium  values.  In  our  model  runs,  we  used  a  grid  spacing  of  3  .11 
correspond! ng  to  the  nominal  resolution  of  existing  airborne  SAR  sys~ 
terns.  At  longer  wave'engths,  the  perturbations  caused  by  a  given 
current  are  smaller  than  for  shorter  wavelengths,  but  the  relaxation 
effects  are  also  reduced.  Th<s  is  illustrated  in  Figure  13  for  the 
ICC  cm  wavelength  case.  The  perturbation  of  these  waves  is  smaller 
than  the  shorten  waves  for  little  or  no  wind,  but  it  dees  net  change 


significantly  as  wind  effects  are  included.  One  word 
snculc  be  included.  Our  "odel  assumes  a  saturated  Phi 


if  caution 
PS  Spectrum; 


past  radar  and  surface  measurements  have  indicated  that  this  assump¬ 
tion  is  justified  for  X-  and  L-:and  (3.2  and  25  cm)  wavelengths;  it 
has  not  been  demonstrated  for  100  cm  waves.  It  should  a'so  be  noted 
that  the  low  wind  speed  simulations  may  not  represent  a  realistic 
situation  due  to  the  lack  of  background  signal. 


A  second  set  of  simulations  were  generated  for  the  500  kiloten 
supertanker  case  corresponding  to  a  'ocation  five  kilometers  behind 
the  ship  (Case  3).  The  peak  lateral  current  velocity  for  this  case 
is  acout  12  cm/s  as  shown  on  Figure  11.  we  again  adjusted  the 

vcrticity  in  our  node’  to  produce  similar  currents.  The  remaining 
parameters  were  fixed  at  values  predicted  by  theory.  The  surface 
current  profile  calculated  for  this  case  '  s  shewn  plotted 
figure  1A.  The  maximum  spectral  perturbations  caused  by  this  current 
field  is  shown  in  Figure  15  as  a  function  of  surface  wavelengtn  and 
wind  speed.  Due  to  the  large  currents  assoc’ated  with  this  case, 
our  model  cculd  not  be  used  for  the  zero  wind  case  or  for  the  two 
snorter  wavelength  examoles.  However,  the  cases  which  were  run 

illustrate  tne  important  points. 

Similar  to  the  Case  1  stucy  presented  above,  the  spectral  oer- 
turoations  caused  by  the  vortex-re  1 ated  currents  is  a  maximum  at 
shorter  wavelengths  for  ve^y  'cw  wind  speeds  and  at  longer  *ave- 
lengt.ns  for  nigner  wind  speeds.  In  contrast  to  tne  Case  1  results 
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MAXIMUM  SPECTRAL  PERTURBATION 
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gure  i5.  Perturbations  Caused  by  tti®  Vortex-Produced  Surface 
Currents  in  Figure  14  as  a  Function  of  Wavelength 
and  Wind  Speed 
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presented  in  Figure  13,  the  negative  perturbation  of  100  cm  waves  is 
very  hign  for  low  wind  speeds  but  decreases  fairly  rapidly  as  the 
wind  increases.  Recall  that  for  Case  1,  the  spectral  perturbation  of 
the  100  cm  waves  remained  nearly  constant  for  the  wind  speeds  con¬ 
sidered.  It  is  also  interesting  to  note  that  the  perturbation  of 
the  25  cm  waves  for  this  case  is  greater  than  tne  Case  1  perturba¬ 
tions  at  1.5  and  3  m/s  wind  speeds,  but  less  than  the  Case  1  pertur¬ 
bations  for  5  and  12  m/s  winds.  These  results  can  be  explained  with¬ 
in  the  context  of  Our  model. 

A:  very  lew  wind  speeds,  the  spectral  perturbation  of  surface 
waves  is  proportional  to  the  surface  current  velocity.  At  very  hign 
wind  speeds,  the  spectral  perturbation  is  proportional  to  the  strain 
rate(du/dy).  which  is  the  spatial  derivative  of  the  surface  current 
profile.  At  intermediate  wind  speeds,  the  spectral  perturbation  is 
proportional  to  some  combination  of  current  and  strain  rate  effects. 
The  wind  speed  at  which  the  strain  rate  becomes  important  depends  on 
the  surface  wavelength  being  considered  and  occurs  at  a  lower  wind 
speed  for  short  waves.  Ih  the  observations  above,  the  surface  cur¬ 

rents  in  Figure  12  were  only  about  one-third  as  large  as  those  in 
-  igure  1*,  but  the  strain  rates  were  about  twice  as  large,  'he  *act 
that  the  perturbation  of  tne  100  cm  waves  in  the  Case  1  simulation 
remain  nearly  constant  indicates  that  the  relative  magnitude  of  the 
current  and  strain  rate  effects  are  approximately  equal.  The 
decrease  in  spectral  perturbation  of  the  100  cm  waves  for  Case  3  as 


the  wind  speeds  increase  indicates  that  the  strain  rate  contribution 


is  becoming  important.  Finally,  the  fact  that  the  perturbation  of 
the  25  cm  waves  is  greater  in  the  Case  1  study  for  wind  speeds  of  5 
and  12  m/s,  even  though  the  surface  currents  were  smaller,  clear’y 


indicates  that  *or  this  wavelength  ana  these  wind  speeds,  the  str 
rate  is  the  dominant  factor. 

A  third  and  final  set  of  simulations  were  generated  for  the 
< i lotcn  supertanxer  case  'Case  3)  corresponding  to  a  location 
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<i1ometar$  behind  the  ship.  The  peak  lateral  current  velocity  fp r 
this  case  is  about  4,5  cm/ s  as  shown  on  Figure  11,  The  orticity  in 
Our  mode  1  was  adjusted  to  natch  these  currents.  The  remaining  para¬ 
meters  were  consistent  with  the  previous  set  of  simulations.  The 
surface  current  profile  for  this  case  is  shown  plotted  in  Figure  16. 
The  spectral  perturbations  caused  by  these  current  fields  is  shown  in 
Figure  17  as  a  function  of  wavelength  and  wind  speed.  These  results 
are  what  we  would  expect  based  on  the  discussion  above.  Both  the 

surface  currents  and  strain  rates  are  only  about  one-third  of  their 
values  in  tne  previous  set  of  simulations.  This  reduces  the  spectral 
pert jrbations  at  octh  low  and  high  wind  speeds. 

A  final  plot  was  produced  which  illustrates  the  spectral  pertur¬ 
bations  caused  by  vortex-related  wakes 'as  a  function  of  distance  aft 
of  the  generating  ship.  This  plot  is  shown  as  Figure  18  and  was 
generated  from  the  Case  3  simulations  discussed  above.  These  pertur¬ 
bations  correspond  to  locations  5  and  10  km  behind  a  supertanker  and 
are  plotted  for  25  and  ICO  cm  wavelengths  and  wind  speeds  ranging 
from  1.5  to  12  m/$.  These  results  indicate  that  the  spectral  pertur¬ 
bations  significantly  decrease  due  to  the  decay  of  the  vortex¬ 
generated  surface  currents.  These  results  indicate  that  25  cm  sur¬ 
face  waves  would  net  oe  detectable  at  either  location  under  12  -'s 
winds.  At  10  km  behind  the  ship,  the  25  cm  waves  would  net  be 
detectable  with  6  or  12  m/s  winds  and  only  marginally  detectable  ‘‘or 
3  m/s.  The  ICO  cm  surface  waves  would  only  be  marginally  detectec  at 
the  10  km  position  for  12  m/s  winds. 

The  simulations  presented  above  have  indicated  several  things. 
At  very  low  wind  speeds,  the  spectral  perturbations  associated  th 
sn i p-generateo  vortices  are  a  maximum  for  shorter  wavelengths,  -s 
wind  speed  increases,  wave  growth  or  relaxation  effects  cause  tne 
■perturbation  of  the  snorter  waves  to  decrease  sharply.  Although  they 
are  net  as  strongly  perturbed  for  a  given  current  pattern  at  ’cw  w-r -d 
speeds ;  longer  surface  waves  are  more  insensitive  to  -el axat  *  on 


Figure  16.  Modelled  Vortex-Produced  Surface  Currents  ler* 
Kilometers  Behind  a  500  Ktloton  Ship  (Case  3} 
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figure  17.  Hovelled  Perturtati  ons  of  the  Vortex-Produced 
Surface  Currents  in  Figure  16  as  a  Function  of 
Wavelength  and  Wind  Speed 
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Figure  18.  Modelled  Perturbations  of  Vortex-Produced  Wakes 

for  Case  3  as  a  Function  of  Wavelength,  Wind  Speed1 
and  Distance  Aft  of  the  Ship 
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effects  and  offer*  the  possibility  of  detection  at  higher  wind  speeds. 
Ideally,  one  would  want  a  sensor  which  would  be  sensitive  to  the 
short  waves  under  low  winds,  and  longer  waves  for  higher  wind  speeds. 
One  possibility  would  be  a  multiple  frequency  SAR  system. 

The  results  of  our  simulations  also  indicate  that  the  type  of 
ship  generating  the  vorticies  also  determines  their  detectability. 
At  low  wind  speeds,  the  spectral  perturbation  of  surface  waves  is 
proportional  to  surface  currents,  which  in  our  model  are  proportional 
to  ship  velocity  and  beam  and  inversely  proportional  to  the  square 
root  of  the  length.  At  higher  wind  speeds,  the  spectral  perturbation 
is  related  to  the  strain  rate.  The  strain  rates  would  be  highest  for 
those  cases  where  the  vorticies  shed  by  the  Ship  have  not  spread 
apart  significantly.  This  would  be  true  for  narrow  ships  with  high 
velocities. 

The  reader  should  be  cautioned  that  the  model  we  have  used  above 
is  only  a  first-order  approach  to  tne  problem.  The  effects  of  the 
rearward  axial  current  should  be  included  as  part  of  extending  the 
model  to  two  dimensions.  We  would  also  like  to  include  SAR  imaging 
considerations  into  the  model  such  as  motion  effects,  scatterer 
coherence  time,  and  variable  resolution. 

2.3  ADDITIONAL  IMAGERY 

Presented  below  are  several  SAR  images  of  ships  and  what  appear 
to  be  turbulent  wake  features.  Very  little  is  known  about  the  ship 
parameters  in  these  images  or  the  surface  wind  conditions  during  the 
data  collections.  This  limits  our  discussions  of  these  images  to  a 
general  qualitative  nature. 

Figure  19  represents  simultaneously  collected  L-  and  X-band  SAR 
imagery  af  a  ship  towing  a  barge.  In  each  image,  a  turbulent  wake 
feature  is  observed  aft  of  the  barge.  In  the  L-band  image,  the  wake 
feature  is  observed  to  have  a  bright  edge  and  a  band  of  reduced 
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Figure  19.  L-  and  X-Band  Optical  SAR  Imagery  of  a  Turbulent 
Wake  Feature  Collected  2  November  1978 
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backscatter  which  gradually  returns  bo  the  background  intensity  about 
5  km  aft  of  tne  barge,  In  tne  X-band  i mag€ *  a  narrow  bright  arm  <s 
observed  and  also  a  very  narrow  da r<  return  area  extending  to  about 
<1  km  behind  the  snip.  We  believe  that  the  wind  direction  for  these 
data  was  towards  the  top  of  the  page.  This  is  bated  on  the  bright 
arms  observed  in  the  imagery  where  we  feel  the  vortex-induced  cur¬ 
rents  and  wind-generated  capillary  waves  are  in  opposition.  We  mace 
a  rough  attempt  at  modelling  these  features  by  assuming  a  current 
pattern  cf  the  same  ^orm  as  Case  1  in  the  previous  section.  In  cu,‘ 
model  -uns,  we  assumed  a  wind  speed  of  only  1  m/s.  The  surface  cur¬ 
rent  profile  anc  spectral  perturbati cns  predicted  by  our  model  are 
shown  in  Figure  20  *o r  the  L-ba*.d  Bragg  waves.  The  modelled  spectral 
perturbations  qualitatively  agree  with  the  image  observation,  that 
is,  a  narrow  line  of  increasec  backscatter  on  the  assumed  upwind  side 
of  the  wake  and  a  broader  region  of  decreased  backscatter  downw.nd. 
The  modelled  perturbations  for  tne  X-band  3ragj  waves  are  presented 
in  Figure  21.  rhese  results  do  not  agree  as  well  with  the  imagery  as 
the  l-oand,  but  nevertheless,  several  points  can  be  made.  The  model 
oces  net  predict  tne  bright  narrow  arm  on  the  assured  upw'nd  s;ce  o* 
the  wake  which  is  observed  in  the  data,  we  do  not  f,jl ly  understand 
tn  is  but  think  that  the  SAS  may  be  sensing  snort  waves  created  by  tne 
break ing  of  longer  waves.  The  model  does  predict  the  narrow  line  c# 
decreased  bacxscatter  observed  in  the  imagery.  rhe  structure  of  tne 
modelled  perturbations  in  Figure  21  indicates  that  there  may  be  cases 
where  two  lines  of  decreased  backscatter  extend  benind  a  snip  with 
the  water  between  them  having  re': axed  back  to  equilibrium.  This  has 
in  fact  been  observed  in  .numerous  cases  and  is  referred  to  as  turbu¬ 
lent  waxe  splitting.  An  example  is  shown  -in  Figure  22. 
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Figure  20.  Estimated  Surface  Current  Profile  and  Modelled 
Spectral  Perturbations  for  L-8and  Bragg  Waves 
Corresponding  to  Figure  19 
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Figure  21.  Estimated  Surface  Current  Profile  and  Modelled 
Spectral  Perturbations  for  X-Band  Bragg  Waves 
Corresponding  to  Figure  19 


Figure  22.  X-Band  Optical  SAR  Image  Indicating  Turbulent  Wake 
Splitting  (Pass  2,  3  November  1978) 
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Figure  23.  Seasat  5AR  Image  of  Suspected  Oil  Tanker  in  the 
Gulf  of  Alaska  (after  Fu  and  Holt,  1982) 
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were  measured  during  this  data  collection.  Using  the  Case  3  current 
profile  corresponding  to  10  km  behind  the  ship,  a  peak  spectral  per¬ 
turbation  of  0.06  is  found  for  these  wind  conditions.  Earlier,  we 
qualitatively  determined  that  a  spectral  perturbation  of  0.15  *  0.05 
was  necessary  for  detection  of  a  wake  feature  in  $AR  imagery.  Recall 
that  Seasat  SAR  data  was  typically  processed  to  four  looks;  this 
reduces  our  detectability  requirements  to  about  0.075  *  0.025  so  our 
model  results  appear  consistent  with  the  image  observations. 

Two  X-oand  SAR  images  of  ships  and  their  wakes  are  presented  in 
Figure  24.  These  images  are  both  from  the  same  pass  of  data  sepa¬ 
rated  by  about  20  km.  The  wind  speeds  during  the  collection  of  these 
data  was  very  low,  in  fact,  visual  observations  during  these  flights 
indicated  that  the  ships  were  the  only  source  of  surface  roughness. 
The  ship  in  Figure  24(a)  is  displaced  from  its  wake  due  to  its  range 
velocity  component.  The  speed  of  this  ship  has  been  estimated  from 
this  displacement  at  about  10  m/s.  The  radar  return  behind  this  ship 
is  concentrated  into  two  narrow  lines  extending  aft  of  the  ship.  In 
Figure  24(b),  the  --adar  return  behind  the  ships  consists  of  a  region 
of  slightly  increased  backscatter.  The  Doppler  displacement  of  these 
ships  indicated  speeds  of  about  6  m/s.  Recall  from  Section  2.1  that 
the  magnitude  of  surface  currents  associated  with  ship-generated 
vortices  is  directly  proportional  to  ship  speed.  We  feel  that  the 
higher  speed  of  the  ship  in  Figure  24(a)  is  creating  higher  surface 
currents  which  concentrate  the  sh ip-qenerated  roughness  along  two 
narrow  lines  corresponding  to  the  vortex  locations. 

Our  final  example  is  an  X-band  optical  image  of  a  large  ship  in 
the  Norwegian  Sea  presented  as  Figure  25.  We  have  estimated  the  wind 
speeds  as  10  m/s  during  the  data  collection.  Observed  on  the  image 
are  the  two  arms  of  the  Kelvin  wake  and  two  bright  narrow  lines 
extending  behind  the  ship.  At  these  wind  speeds,  the  perturbation  of 
X-band  Bragg  waves  due  to  wave-current  interactions  is  close  to  zero. 
One  possibility  is  that  the  wake  features  in  this  image  are  due  to 
the  generation  of  Bragg-scale  waves  by  the  breaking  of  longer  waves. 
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Figure  24.  X-Band  Optical  SAR  Imagery  Collected  During  Very 

Low  Wind  Conditions  of  (a)  a  Relatively  Fast  Moving 
Ship,  and  (b)  Slow  Moving  Ships 
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Figure  25.  X-Band  Optical  SAR  Image  Indicating  Wave  Breaking 
Due  to  Vortex  Current  Interactions 
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2.4  RESOLUTION  EFFECTS 

In  addition  to  the  radar  wavelength  dependence  of  SAR  ship  wake 
detection,  the  resolution  of  the  SAR  system  being  considered  will 
also  determine  what  ship  wake  components  will  be  detected  and  to  what 
degree.  To  study  the  effects  of  changing  resolution,  a  digital 
L-band  ship  wake  image  collected  during  the  1983  Georgia  Strait 
Experiment  was  reprocessed  with  a  reduced  bandwidth  corresponding  to 
Seasat  parameters.  The  original  full  bandwidth  image  is  shown  in 
Figure  26(a).  This  image  has  a  resolution  of  about  4  m  in  both  slant 
range  and  azimuth.  Visible  in  the  image  are  the  turbulent  wake, 
stern  waves,  and  both  arms  of  the  Kelvin  envelope.  For  a  general 
discussion  of  these  various  components  refer  to  Lyden  et  al .  (1985a). 
This  image  was  reprocessed  with  reduced  bandwidths  corresponding  to 
one-look  Seasat  resolutions  of  8  x  7  m  in  slant  range  and  azimuth, 
respectively,  and  is  presented  as  Figure  26(b).  At  these  reduced 
resolutions,  only  the  turbulent  wake  is  visible,  and  its  detect¬ 
ability  has  dropped  relative  tG  the  full  resolution  image. 

A  comparison  of  actual  Seasat  imagery  with  simultaneously- 
collecteo  aircraft  data  of  a  ship  wake  is  shown  in  Figure  27,  These 
data  were  collected  in  the  Georgia  Strait  in  1978.  Stern  waves  are 
faintly  visible  in  the  aircraft  data  as  well  as  a  clearly  defined 
turbulent  wake.  Only  a  very  faint  turbulent  wake  is  visible  in  the 
Uaasat  data.  A  comparison  of  these  data  with  the  aircraft  and 
aircraft-simulated  Seasat  data  in  Figure  26  reveals  several  things. 
In  both  cases,  stern  waves  were  visible  in  the  aircraft  data  but  not 
in  the  Seasat  case.  Turbulent  wakes  are  also  clearly  visible  in  the 
aircraft  data,  but  are  not  as  clearly  visible  in  the  Seasat  or 
Seasat-simulat.ed  data.  This  comparison  clearly  illustrates  the 
effects  of  resolution  on  the  detectability  of  subtle  patterns,  and 
also  shows  that  these  effects  can  be  simulated  by  appropriate  proces¬ 
sing  of  high  resolution  aircraft  data. 
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(a)  L-B AND  DIGITAL  AIRCRAFT  DATA 


(b)  SEASAT-SIMULATED  DATA 


Figure  26.  t-Band  Digital  Aircraft  SAR  Data  Processed  to 
(a)  Full  Resolution,  and  (b)  Simulated  Seasat 
Resol uti ons 
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The  ability  to  simulate  satellite  ship  wake  SAR  imagery  using 
existing  aircraft  data  is  important  for  several  reasons.  Simulations 
of  this  type  can  be  used  to  determine  resolution  requirements  for  the 
detection  of  various  ship  wake  components.  An  example,  would  be  to 
determine  the  resolution  necessary  to  detect  the  stern  (transverse) 
waves  visible  in  Figure  26(a).  This  could  be  accomplished  by 
repeatedly  processing  this  data  at  coarser  resolutions  until  trie 
stern  waves  disappear.  These  simulations  can  also  be  used  to  predict 
the  ship  wake  signatures  which  can  be  expected  from  future  satellite 
systems. 

2.5  DETECTION  APPLICATIONS 

The  final  study  we  performed  as  part  of  our  ship  wake  analysis 
was  to  examine  the  Fourier  transform  of  a  ship  wake  image  for  unique 
characteristics  which  could  then  be  used  in  an  automated  detection 
algorithm.  The  ship  wake  data  we  used  in  this  study  was  collected 
during  the  1983  Georgia  Strait  Experiment  and  is  presented  in 
Figure  28.  The  wake  feature  in  this  image  is  composed  of  two  bright 
narrow  arms  which  are  modulated  by  the  stern  waves.  This  image  was 
fast  Fourier  transformed  and  the  resulting  spectrum  is  displayed  in 
Figure  29.  These  results  are  displayed  in  wavenumber  space  so  that 
the  center  of  this  image  corresponds  to  zero  frequency  or  DC,  with 
increasing  distance  from  this  point  representing  shorter  wavelengths. 
The  horizontal  and  vertical  scales  on  Figure  29  are  different  because 
we  transformed  a  rectangular  rather  than  a  square  region.  Several 
features  are  visible  in  the  transform.  Two  lines  of  low-frequency 
energy  are  present  corresponding  to  and  aligned  perpendicular  to  the 
narrow  arms  in  the  ship  wake  image.  These  are  what  would  be  expected 
if  the  arms  in  the  image  could  be  modelled  as  solid  lines.  The  spec¬ 
tral  resolution  of  this  transform  allows  these  arms  to  be  dif¬ 
ferentiated.  There  is  also  a  well-defined  spectral  component  which 
is  nearly  vertical ly-al ignea  with  a  wavelength  of  about  35  m.  This 
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Figure  29.  Fast  Fourier  Transform  Magnitude  Spectrum  of  SAR 
Image  Shown  in  Figure  28 
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is  the  mcd u  1  at  1  on  we  have  associated  with  the  3$  m  wavelength  stern 
waves,  These  waves  travel  at  the  same  velocity  as  the  ship  generat¬ 
ing  them,  thus,  they  are  very  mono chromatic  fur  a  constant  velocity 
ship.  With  some  pre-filtering,  to  remove  noise  effects,  this  com¬ 
ponent  could  be  enhanced.  Since  monochromatic  waves  at  "engths  of 
30-60  m  are  not  expected  in  the  open  ocean,  this  component  being 
present  in  Fourier  transformed  SAR  data  may  indicate  sh io-generated 
stern  waves. 

A  subset  of  the  ship  wake  image  discussed  above  is  shown  in 
Figure  30.  This  subset  was  extracted  from  the  image  shown  as 
Figure  27  and  its  Fast  Fourier  Transform  (FFT)  was  calculated.  The 
spectrum  of  this  transform  is  presented  in  Figure  31,  Several 
interesting  features  are  present  in  this  spectrum.  The  wake  arms  are 
not  resolved  as  they  were  in  Figure  29.  This  is  because  spectral 
resolution  is  deg'aded  as  the  area  being  transformed  Is  decreased. 
Again,  there  is  a  significant  component  at  about  40  m  wavelength 
which  is  associated  with  the  stern  waves. 

The  work  described  above  is  very  preliminary  in  nature. 
Additional  work  in  this  area  would  include  modelling  the  expected 
transforms  of  ship  wake  features,  and  seeing  how  the  modelled  results 
compare  with  transforms  of  ratural  phenomena.  This  would  provide 
insight  into  the  uniqueness  of  the  ship  wake  features.  Other  tech¬ 
niques  besides  the  Fourier  transform  should  also  be  investigated  for 
automated  detection  applications. 
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Figure  30.  Subset  of  Ship  Wake  Image  Stown  in  Figure  28  Used 
in  Detection  Studies 
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Figure  31.  Fast  Fourier  Transform  Magnitude  Spectrun  of  SAR 
Image  Snown  in  Figure  30 
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CONCLUSIONS  and  recommendations 


The  primary  purpose  of  this  investigation  was  to  study  the  darn 
wakes  observed  in  SAR  imagery  behind  moving  ships,  A  first  attempt 
has  been  made  at  using  a  hydrodynamic  model  *or  surface  currents 
associated  with  ship-generated  vortices  as  input  into  a  wave-current 
interaction  model  which  predicts  a  spectral  perturbation  of  the  wave¬ 
length  under  consideration.  These  simulations  were  performed  at 
wavelengths  corresponding  to  the  Bragg  waves  for  various  radars. 
Spectral  perturoations  of  their  Bragg  waves  are  what  radars  are  sen¬ 
sitive  to.  For  the  cases  considered,  the  simulated  wake  signatures 
generally  agreed  to  these  observed  in  the  actual  data.  Our  results 
indicate  several  things. 

The  Spectral  perturbations  associated  with  a  vortex-rel ated  sur¬ 
face  current  pattern  are  larger  at  shorter  wavelengths  than  at  longer 
wavelengtns.  Therefore,  at  low  wind  speeds  where  wind  effects  are 
negligible,  a  vortex-related  surface  current  pattern  would  be  most 
detectable  by  a  radar  wh'cn  was  sensitive  to  shorter  waves.  In 
higher  wind  conditions,  the  shorter  wavelengths  regenerate  mere 
quickly  than  the  longer  wavelengths.  Therefore,  although  the 
vortex-rei ated  surface  currents  have  a  greater  effect  on  the  snorter 
waves,  these  waves  quickly  return  to  saturation.  This  return  to 
equilibrium  occurs  over  a  faction  of  a  typical  SAR  resolution  cell 
thereby  masking  any  current-related  features  in  the  imagery. 
Although  longer  surface  waves  are  -not  as  sensitive  to  surface  cur¬ 
rents,  they  regenerate  more  slowly  in  response  to  wind.  ”herefore, 
at  higher  wind  speeds,  a  vortex-re  1 ated  surface  current  pattern  wcu'd 
be  most  detectab’e  by  a  radar  which  was  sensitive  to  longer  waves. 
’deal1;/,  ere  would  want  a  radar  which  wou’d  be  sensitive  to  the  short 
waves  under  lew  winds,  and  longer  waves  for  higher  wind  speeds .  Ore 
possio':l;ty  would  be  a  multiple  frequency  SAR  system. 
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The  results  of  our*  simulations  also  indicate  that  the  type  of 
ship  generating  the  vortices  determines  their  detectao i 1 ity.  At  lew 
wind  speeds,  the  spectral  perturbation  of  surface  waves  is  propor¬ 
tional  to  surface  currents,  which  in  our  model  are  proportional  to 
ship  velocity  and  beam  and  inversely  proportional  to  the  square  root 
of  the  length.  At  higher  wind  speeds,  the  spectral  perturbation  is 
related  to  the  strain  rate  which  is  the  spatial  derivative  of  the 
Surface  current  profile.  Strain  rates  would  be  hignest  for  those 
cases  where  the  vortices  shad  by  the  sh'p  nave  not  spread  apart  sig¬ 
nificantly.  This  would  be  true  for  narrow  ships  with  hi  gn 
velocities . 

Recommendations  for  future  work  in  this  area  would  be  to  include 
the  rearward  axial  current  into  our  model  and  to  extend  this  model  to 
perform  two-dimensional  simulations.  This  rearward  current  is 
initially  smaller  than  the  vortex-rel ated  lateral  currents,  but 
decays  mere  slowly  with  time.  We  would  also  like  to  include  the 
effects  of  coherence  time,  scatterer  motion,  and  variable  resolution 
in  our  model . 

Other  analyses  wnich  were  performed  as  part  of  this  study  ;n- 
cljded  examining  the  **e '  at  ipnsh '  p  between  SAR  resolution  anc  toe 
v'sipility  of  various  ship  wake  components.  This  was  done  by  proces¬ 
sing  aircraft  SAR  data  of  a  ship  waxe  with  a  reduced  bandw’dtn  cor¬ 
responding  to  the  Seasat  SAR.  A  comparison  of  the  two  processings 
indicated  that  many  of  the  srrp  wake  components  were  lost  at  the 
coarser  resolutions.  A  comparison  of  simu’ taneously  collected  air¬ 
craft  ana  Seasat  SAR  aata  indicated  that  our  Seasat  simulation  pro¬ 
duced  reasonable  results.  Js^ng  existing  aircraft  data,  we  can  ^se 
tnis  approacn  to  predict  wnat  ship  wake  features  future  satellite 
SARs  snould  oe  ab’e  to  detect  as  well  as  help  to  define  resc'ut-cr 
requirements  for  future  SAR  design. 
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The  final  pant  c f  this  study  examined  the  Fourier  transform  of 
$AR  ship  m ake  imagery  for  unique  information  'which  could  he  used  in 
automated  detection  algorithms.  Results  of  this  study  were  incon¬ 
clusive  hut  indicated  that  there  is  some  information  in  the  ship  wake 
image  transform  which  may  be  unique  and  be  able  to  be  used  for 
detection  purposes.  Recommendations  for  future  research  include 
model 'ing  observed  ship  wake  signatures  and  their  transforms  and  com¬ 
paring  these  to  targets  which  would  be  expected  to  occur  naturally. 
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APPENDIX  A 

SAR  DETECTABILITY  OP  SURFACE  WAVE  PERTURBATIONS 

The  detectability  of  vortex  wake-induced  current  patterns  by  SAR 
depends  not  only  on  the  spectral  perturbation  of  the  Bragg  wave¬ 
lengths  caused  by  the  currents,  but  also  on  the  site  of  the  perturbed 
region  and  the  resolution  and  speckle  characteristics  of  the  SAR 
system.  Presented  below  is  a  brief  discussion  on  how  these  various 
factors  influence  the  detectability  of  vortex-related  patterns  in  SAR 
imagery.  Several  image  simulations  are  presented  which  help  illu¬ 
strate  these  effects. 

The  model  we  have  employed  in  our  vortex  wake  simulations  calcu¬ 
lates  the  spectral  perturbation  at  a  specified  wavelength  due  to 
interactions  with  vortex-induced  surface  currents  (see  Section  2.1). 
Spectral  perturbations  of  the  Bragg  wavelengths  are  assumed  to  be 
directly  proportional  to  changes  in  radar  backscatter.  The  detect¬ 
ability  of  these  backscatter  changes  in  a  SAR  image  is  dependent  on 
the  magnitude  of  the  spectral  perturbation,  the  extent  of  the  per¬ 
turbed  region  relative  to  the  SAR  resolution,  and  the  speckle  char¬ 
acteristics  of  the  SAR  system. 

Simulated  SAR  images  were  produced  for  the  spectral  perturbation 
pattern  shown  in  Figure  A-l.  This  pattern  consists  of  different 
width  lines  with  spectral  perturbations  ranging  from  0.0  on  the 
extreme  left  of  the  pattern  to  0.5  on  the  extreme  right  side.  The 
vertical  width  of  these  features  ranges  from  two  pixels  at  the  top 
of  the  pattern  to  16  pixels  at  the  bottom.  Intuitively,  the  detect¬ 
ability  of  these  features  should  increase  with  size  and  spectral 
perturbation. 


The  simulated  SAR  image  presented  in  Figure  A-2  was  generated 
assuming  a  resolution  equal  to  2  x  2  pixels  and  one-look  speckle 
statistics.  Recall  that  speckle  is  a  multiplicative  noise  effect 
which  arises  due  to  the  coherent  nature  of  SAR  (Porcello,  et  al.. 
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Figure  A-2.  One-Look  Simulated  SAR  Image  Illustrating  the 

Detectability  of  Wake  Features  as  a  Function  if 
Spectral  Perturbation  and  Width 
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1976).  This  image  indicates  sever'd!  things.  The  narrowest  perturbed 
region  is  not  visible  in  the  image,  even  for  the  highest  perturbation 
values.  The  wider  features  are  detectable  at  various  perturbation 
values  with  the  widest  features  being  visible  at  the  smallest  pertur¬ 
bations.  To  accurately  determine  at  what  perturbation  a  feature  in 
this  image  becomes  detectable,  the  image  should  be  examined  from  left 
to  right  with  the  right  hand  portion  (Urge  perturbation  region) 
covered.  From  Figure  A-2,  it  appears  that  a  perturbation  of  about 
0.15  *  0.05  is  detectable  for  line  widths  greater  than  about  6  reso¬ 
lution  cells.  A  four-look  simulation  of  the  pattern  in  Figure  A-l  is 
shown  in  Figure  A-3.  The  resolution  for  this  simulation  was  again 
assumed  to  be  2  x  2  pixels.  The  standard  deviation  of  the  speckle 
in  the  four-look  data  is  reduced  by  a  .factor  of  two  over  the  one-look 
data.  This  improves  the  detectability  of  the  various  features.  The 
narrowest  feature  is  now  visible  at  high  perturbations.  The  minimum 
detectable  perturbation  now  appears  to  be  about  0.075  *  0.025  for 
image  features  wider  than  six  resolution  cells. 

Models  of  detectability  based  on  matched  filter  theory 
(Lyzenga,  1935)  suggest  that  the  index 


dx  dy 


1/2 


!A-1) 


where  (ao/o  )  represents  the  pattern  strength  and  o  o  represent 

y 

the  SAR  resolution,  can  be  used  to  estimate  detectability. 
Specifically,  a  value  of  D  in  the  order  of  15  has  postulated  as  a 
threshold  for  visual  detectability  of  internal  wave  patterns  in  cases 
where  speckle  noise  dominates  the  background  visibility. 


For  a  linearly  increasing  pattern  strength  as  used  in  the  simula¬ 
tions  presented  here,  the  detectability  index  for  the  portion  of  the 
pattern  to  the  left  of  the  position  X  may  be  written  as 


0(X) 


(A-2) 
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Figure  A-3.  Four-Look  Simulated  SAR  Image  Illustrating  the 
Detectability  of  Wake  Features  as  a  Function  of 
Spectral  Perturbation  and  Width 
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where  W  is  the  pattern  width,  &  is  the  SAR  resolution,  and  X  is 

0 

the  position  at  which  tala  *  1.  Comparison  of  this  index  with  the 
simulated  images  shows  that  the  patterns  become  visible  where  0^4, 
which  is  the  theoretical  limit  for  detectability  using  a  matched 
;  ilter.  The  reason  for  the  lower  apparent  threshold  of  detectability 
for  these  patterns  (as  opposed  to  internal  wave  patterns)  is  their 
simpler  shape.  That  is,  a  pattern  consisting  of  a  single  line  is 
more  easily  detectable  than  a  more  complicated  pattern  having  the 
same  perturbation  strength  and  area. 

Our  final  simulation  also  corresponded  to  the  perturbation  pat¬ 
tern  shown  in  Figure  A-l,  b'.’t  was  performed  assuming  a  resolution 
equal  to  one  pixel  and  one-look  speckle  statistics.  The.  ©suiting 
image  is  shown  in  Figure  A-4.  This  image  appears  very  similar  to  the 
four-look  simulation  in  Figure  A-3,  even  though  the  variance  is  twice 
as  high.  This  is  a  result  of  the  viewer's  eye  performing  the  non¬ 
coherent  integration  (multi-looking)  of  the  imagery.  The  image  shown 
in  Figure  A-4  would  be  identical  to  Figure  A-3  if  it  was  smoothed 
using  a  2  by  2  pixel  boxcar  filter.  In  other  words,  this  image  is 
the  cne-1nok  equivalent  to  Figure  A-3  (i.e.,  both  images  have  the 
same  bandwidth).  The  four-look  eauivalent  to  Figure  A-2  would  have  a 
resolution  of  4  by  4  pixels. 

In  summary,  we  have  performed  the  simulations  presented  above  to 
show  the  effects  of  target  size,  SAR  resolution,  and  speckle,  on  the 
detectability  of  features  with  varying  spectral  perturbations.  These 
simulations  show  that  features  become  more  detectable  with  increasing 
size,  increasing  spectral  oerturbation,  and  finer  resolution  or  in¬ 
creasing  number  of  looks. 
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Figure  A-4.  Higher  Resolution  One-Look  Simulated  SAR  Image 

Illustrating  the  Detectability  of  Wake  Features  as  a 
Function  of  Spectral  Perturbation  and  Width 
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